Intermittent implosion and pattern formation of trapped Bose-Einstein condensates 

with attractive interaction 



Hiroki Saito and Masahito Ueda 
Department of Physics, Tokyo Institute of Technology, Tokyo 152-8551, Japan 
and CREST, Japan Science and Technology Corporation (JST), Saitama 332-0012, Japan 

(February 1, 2008) 



O 
O 
O 

> 
O 

(N 

> 
cn 
o\ 
cn 
(N 
O 
O 

o 



I 

o 
o 



The collapsing dynamics of a trapped Bose-Einstein con- 
densate (BEC) with attractive interaction are revealed to ex- 
hibit two previously unknown phenomena. During the col- 
lapse, BEC undergoes a series of rapid implosions that oc- 
cur intermittently within a very small region. When the sign 
of the interaction is suddenly switched from repulsive to at- 
tractive, e.g., by the Feshbach resonance, density fluctuations 
grow to form various patterns such as a shell structure. 
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Bose-Einstein condensation (BEC) of trapped atomic 
vapor has been realized in ^''Rb 0], ^■^Na and 
^Li iQ . The last species is unique in that it has a negative 
s-wave scattering length, implying that the interactions 
between atoms are predominantly attractive. It has been 
believed that in a spatially uniform system, such atomic 
vapor would collapse into a denser phase. However, when 
the system is spatially confined and when the number 
of BEC atoms is below a certain critical value iVc, the 
zero-point motion of the atoms serves as a kinetic ob- 
stacle against collapse, allowing a metastable BEC to be 
formed Just below N^, BEC may collapse via 

macroscopic quantum tunneling ]T^- [T8| , p^ , p^Jl3| , and 
above N^, it is predicted that the collapse will oc- 
cur not globally, but only locally near the center of BEC 
where the atomic density exceeds a certain critical value. 
Inelastic collisions also lead to the decay of BEC in re- 
gions where the atomic density is very high |^]. 

In current experiments Q, there are abundant above- 
condensate atoms that replenish the lost atoms, allowing 
BEC to grow again. We may therefore expect coUapse- 
and-growth cycles of BEC to occur. Various mecha- 
nisms that give rise to these oscillations have been dis- 
cussed [^0|- p2|JT^ ] . And indeed, recent experiments ||2^ ] 
have suggested the occurrence of dynamic collapse-and- 
growth cycles of BEC, but the results have neither fa- 
vored nor excluded any one of these possible mechanisms. 

The s-wave scattering length of atoms can be varied 
using the Feshbach resonance suggesting that not 
only the strength of the interaction but also its sign can 
be controlled. Using this technique, the sign of the inter- 
action has recently been successfully switched from pos- 
itive to negative in BEC of ^^^Rb Kagan et al. ^] 
have discussed a global collapse of the condensate in such 
situations. 



In this Letter, we predict two new phenomena associ- 
ated with the collapse of BEC. One is intermittent implo- 
sion, in which the local collapses occur in rapid sequence 
at the center of the condensate. This phenomenon is 
caused by competition between the attraction of atoms 
towards the trap center and the loss of atoms by inelas- 
tic collisions. While our analysis is based on the theory 
developed by Kagan et al. [|2^, this phenomenon is dif- 
ferent from the coUapse-and-growth cycles and other fine 
structures predicted by them. The other prediction is 
that of pattern formation in the atomic density, follow- 
ing a sudden switch in sign of the interaction from re- 
pulsive to attractive. This phenomenon occurs because 
density fluctuations caused by the change in sign of the 
interaction grow and self-focus due to the attractive in- 
teractions. 

We consider a system of Bose-condensed atoms with 
mass m and s-wave scattering length a, confined in 
a parabolic potential. The transition amplitude of 
the system from the initial state ipi to the final 
one tpf is expressed in terms of path integrals as 
J^f 2:>VPV*e^^['^''^'l, where S[tp,^*] is given by 



S[ip, V'*] = Noh dr dt 



(1) 



Here the length, time, and ■0 are normalized in units 
of c?o = {h/mujoy^^ , ujq^, and {No/dl)^/'^, respectively, 
and g = AnNoa/do is the dimensionless strength of the 
interaction. The wave function is then normalized to 
unity. The most probable Feynman path that makes the 
action (|l|) extremal satisfies the Gross-Pitaevskii (CP) 
equation |27[ 
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(2) 



The metastability of BEC with attractive interac- 
tions may be understood by the Gaussian approxima- 
tion If we approximate the wave function as 
having a Gaussian form, the size of BEC — R — obeys 
the equation of motion R — —dVes/dR, where the ef- 
fective potential has the form VcS = 3(i?~^ -I- i?^)/2 -|- 
7i?-3/2 with 7 = The effective potential Fcff has 

a local minimum when I7I < 7c = 8 ■ 5~^/*; therefore, 
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the metastable BEC is formed when there are less than 
the critical number of atoms Nc = -^^-§^7c- We have 
performed a numerical integration of the GP equation 
and have confirmed that the Gaussian approxima- 
tion well describes the dynamics of BEC. However, the 
approximation breaks down when a rapid implosion takes 
place. Since we are interested in the behavior of the local 
implosion, we numerically solve the GP equation without 
resorting to the Gaussian approximation. 

Because the peak density grows very high once the 
collapse begins, we must include in the GP equation 
the atomic loss due to inelastic collisions. Following the 
treatment in Ref . , we employ the GP equation with 
loss processes as 

4t^ = -^VV + y + - ^ ( yl^l' + tI^I') 

(3) 

where L2 and L3 denote the two-body dipolar and three- 
body recombination loss-rate coefficients, respectively. 
The two-body (three-body) loss-rate coefficients must 
be divided by two (six) because of Bose statistics [ p8| . 
We assume that the atoms and molecules produced by 
inelastic collisions escape from the trap without affect- 
ing the condensate. Taking = 27r x 144.5 Hz |2^] 
and the loss-rate coefficients from Refs. p^^, we have 
L2 = 3.7 X lO-'^iVo and L3 = 2.9 x lO-^^iV^ To integrate 
the GP equation, we employ the finite difference method 
with the Crank-Nicholson scheme ||]. Since the implo- 
sion is extremely rapid, we very carefully controlled the 
stepsize to avoid error propagation during the implosion. 

Figure ^ shows the time evolution of the peak height of 
the wave function {ip^r = 0,t)\ (solid curve), the number 
of BEC atoms NQ{t) (dashed curve), and the absolute 
squared overlap of the wave function with the initial one 
|/dn/;*(r,0)V'(r,i)P (dotted curve) for 7Vo(0) = 1260, 
which is slightly (0.7 %) greater than Nc- This gives an 
estimate of the "condensate fraction" that is measured by 
the absorption or phase-contrast imaging (see discussions 
below). We first prepared BEC in a metastable state 
that lay just below the critical point, and then increased 
\a\ (or tightened the trap potential) so that No\a\/do ex- 
ceeded its critical value. In the early stage of the collapse, 
the atomic density increases very slowly and the inelastic 
collisions are unimportant. At t ~ 2.87 a rapid implo- 
sion breaks out, which is blown up in the inset of Fig. ^. 
If the atomic loss were not included, the peak density 
would grow unlimitedly, and the implosion would occur 
only once. With the atomic loss, however, the implosion 
stops in a very short time, and the peak density shows a 
pulse-like behavior. This implosion occurs intermittently 
several times, and with each implosion several tens of 
atoms are lost from the condensate. As a result, the 
number of atoms decreases in a stepwise fashion, even- 
tually reaching approximately 78 % of its initial value 



No. 

Since the coUisional loss of atoms takes place where the 
atomic density is extremely high, the atomic loss primar- 
ily occurs within a very localized spatial region r 0.01. 
The atomic loss is predominantly due to three-body re- 
combination rather than two-body dipolar decay, since 
the loss rate of the former is proportional to the cube 
of the atomic density. In fact, if the two-body loss is 
ignored, intermittent implosions occur (data not shown). 

The duration of the spike shown in Fig. |l| is typically 
At ~ 10~^. The use of mean-field theory with such 
rapid dynamics can be justified as follows. The mean- 
field approximation is applicable for time scales longer 
than 1/gn, which is of order 1 in the metastable state of 
BEC (Note that the time is measured in units of ujq^). 
In the region of implosion, the density n becomes more 
than 10'* times as high as that of the metastable state, 
i.e., l/gn < 10^^, and thus the mean-field theory is valid 
even with the rapid implosion if the relevant time scale 
is longer than 10^"* — which is the case for the situation 
shown in Fig. |l|. The gas parameter na^ is, on the other 
hand, ~ 10~^ in the region of implosion, and the atoms 
are still acting in a weakly interacting regime. 

The pulse-like behavior of the peak atomic density may 
be interpreted as follows. Initially, the condensate has 
a negative pressure due to attractive interactions and 
shrinks towards the central region. When the peak den- 
sity becomes ~ L3|V'|'*/12, i.e., ~ 12g/i3, the 
collisional loss rate of the atoms becomes comparable to 
the accumulation rate of atoms at the center. Since the 
kinetic and interaction energies depend on the atomic 
density and its square, respectively, the total energy in- 
creases upon the loss of atoms . The atoms near the 
center of BEC thus acquire outward momentum, and the 
pressure becomes positive. The change in sign of the 
pressure may be qualitatively explained by the Gaussian 
approximation. Initially —dVcs/dR ~ 3i?~'^ — ||7|i?~''' < 
0, which corresponds to negative pressure; later, however, 
the value becomes positive when the number of atoms (or 
I7I) decreases. After the implosion, inward flow outside 
the region of the implosion replenishes the peak density, 
turning the sign of the pressure again to negative, which 
induces the subsequent implosion. 

When the inward flow is insufficient to reverse the 
sign of the pressure to negative, implosion ceases and 
the atoms are pushed outwards. This phenomenon was 
predicted in Ref. |^2|, and has recently been observed at 
JILA psf l as an atom burst emanating from a remnant 
condensate. According to the estimation in the previous 
paragraph and our numerical analysis [ p^ , the energy 
scale of implosion and subsequent explosion is propor- 
tional to /L3. In the case of ^Li, the mean energy of an 
ejected atom is ~ 80 /iK The atoms and molecules 
are also scattered by three-body recombination, in which 
the release energy is of the order 1 mK. The experimental 
signature of the intermittent implosion should be a series 
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of bursts of atoms and molecules produced by the above 
two mechanisms. 

In the Rice experiments , it has been observed that 
the number of BEC atoms reduces to 10 '--^ 20 % of 
the critical number after the collapse, which is smaller 
than our theoretical evaluation (dashed curve in Fig. |^). 
In Ref. |2^, however, only atoms around the peak of 
the bimodal distribution are collected as the number of 
BEC atoms, while the part of BEC atoms that expands 
broadly following the implosion is hidden in the ther- 
mal cloud. The number of BEC atoms that was mea- 
sured experimentally is roughly estimated from the abso- 
lute squared overlap of the wave function with the initial 
metastable one (dotted curve in Fig. |l|) which decreases 
to below 0.2. Thus our result is consistent with the ob- 
servation of Ref. |2^ . 

The intermittent implosion should be distinguished 
from the two types of oscillations discussed in Ref. [p2| , 
i.e., the coUapse-and-growth cycles and the piecemeal col- 
lapses originating from oscillations of the entire conden- 
sate. The mechanism that causes the intermittent implo- 
sion is quite different from that of those oscillations, as 
discussed above. No intermittent implosion is discussed 
in Ref. ||2|] , as the collisional loss rate used there is much 
larger than that used in this Letter. 

The scenario for the collapse of BEC in the presence 
of condensate growth is, therefore, as follows. The BEC 
grows by being fed by the above-condensate atoms, and 
when Nq exceeds Nc, an implosion occurs that has the in- 
termittent structure shown in Fig. |l|. Some of the atoms 
are lost in the implosion, but they are subsequently re- 
plenished, giving rise to the coUapse-and-growth cycles. 
During these cycles, small collapses occur with period 
^ uJq^ due to oscillations of the entire condensate. These 
small collapses also have intermittent structures, which 
we have confirmed by large-scale numerical simulations. 

We next consider the case in which BEC with repul- 
sive interaction is prepared, and the sign of the inter- 
action is then suddenly switched to attractive. Such a 
situation has recently been realized at JILA ||2^ using 
the Feshbach resonance. Figure || (a) shows the time 
evolution of the wave function, where we assume that at 
t = one million ^^Na atoms having an s-wave scatter- 
ing length of a — 2.75 nm are condensed in the trap with 
frequency ujq = 100 x 2tt s~^. We then change the s-wave 
scattering length to a negative value of —1 nm. We use 
the loss-rate coefficients provided in Ref. pi[, which are 
L2 = 1.7 X lO-^iVo and L3 = 1.1 x 10^N§. When 
the interaction is switched to attractive, the atoms be- 
gin to compress, as the initial wave function has been 
expanded due to the repulsive interaction. The inward 
flow gives rise to a ripple, which then grows up to be a 
series of pulses, as shown in Fig. || (a) (the wave func- 
tion at t = 0.79 is multiplied by 0.1). The growth of the 
density fluctuations can be attributed to the attractive 
interaction, as the atoms tend to accumulate where the 



density is high. The ripple is caused by the momentum 
acquired by the atoms sliding down the trap potential. In 
fact, it can be seen in Fig. ^ (a) that the spaces between 
adjacent pulses near the trap center are smaller than the 
outer spaces. Since the potential energy becomes the ki- 
netic energy as h'^k^/2m ^ muj^R^ where R is the 
initial dimension of the wave function, the wavelength of 
the ripple is estimated to be A/do 2t: da/R. 

The inset in Fig. || (a) shows a gray-scale image 
of the column density integrated along the z axis, 
^'^^\^p{x,y,z)\^dz, sX t = 0.79. This quantity is pro- 
portional to the optical thickness in both absorption and 
phase-contrast imaging, where the laser light propagates 
along the z axis. The concentric circles represent the 
formation of a shell structure in the atomic density. 

The shell-structure formation is due to the instability 
of the initial atomic distribution and the self-focusing ef- 
fect, while the collisional loss is unimportant for that for- 
mation since the atomic density is not so high at the early 
stage of collapse when this phenomenon appears. When 
implosion begins, however, the collisional loss of atoms 
begins to play an important role. Figure |2| (b) shows the 
time evolution of the peak height of the wave function 
\ip{r = 0, f)| (solid curve) and the number of BEC atoms 
(dashed curve). The peak density initially decreases be- 
cause the atoms near the center are attracted towards 
the innermost shell, as shown in Fig. |^ (a). The shells 
move inward, and the first implosion occurs when the in- 
nermost shell arrives at the center of the trap. In Fig. |^ 
(b), five implosions caused by the arrivals of the shells 
at the center of the trap are shown. The number of lost 
atoms associated with such implosions becomes larger for 
the outer shell, since the number of atoms contained in 
each shell is proportional to the square of its original ra- 
dius. The velocity of shells moving inward is roughly 
determined by the free motion of atoms, and the collapse 
time is on the order of uj^^ . Between these implosions, 
the smaller intermittent implosions discussed above oc- 
cur. During each intermittent implosion, several tens of 
atoms are lost, a loss that cannot be discerned in Fig. ^ 
(b) because the total number of atoms is by far greater. 

In the case of an axially symmetric trap, the pattern 
of the atomic density arising from the change of the in- 
teraction is sensitive to the asymmetry of the trap. We 
performed numerical simulations and found that various 
patterns can be formed. For a pancake-shaped trap, the 
atomic motion associated with the change in the inter- 
action is larger in the axial direction than in the radial 
direction. As a result, the ripple arises in an axial direc- 
tion, leading to a layered structure. For a cigar-shaped 
trap, the ripple arises in a radial direction, leading to 
a cylindrical shell structure. These structures undergo 
a complicated evolution and show various patterns such 
as rings and clusters. Intermittent implosions also oc- 
cur for an axially symmetric trap. The details of these 
phenomena will be reported elsewhere. 
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While we have analyzed the specific examples of ^Li 
and ^^Na, the results should be valid for other atomic 
species in which the s-wave scattering length can be var- 
ied, since the parameters g and L3 can be controlled by 
choosing a and Nq. With other values of g and L3, the be- 
haviors may be qualitatively different. When L3 is much 
larger than the value used here, no intermittent implo- 
sion occurs. Formation of a shell structure depends on 
the value of g that is proportional to iVo. For instance, in 
the situation of Fig. ^ when the initial number of atoms 
is A^o = 10^, the number of shells is reduced to two, and 
when Nq — 10*, no shell structure appears. 

In conclusion, we have predicted two new phenomena 
concerning the collapsing dynamics of BEC with attrac- 
tive interactions: intermittent implosions and pattern 
formation in the atomic density. The intermittent im- 
plosion occurs very rapidly compared with the time scale 
of the trap frequency, and in a very localized region com- 
pared with the characteristic size of the trap. When the 
sign of the interaction is suddenly switched from repulsive 
to attractive, the atomic density forms a shell structure 
in a spherically symmetric trap, and various patterns are 
formed for an axially symmetric trap. 
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FIG. 1. Time evolution of the wave function \tp{r = 0,t)\ 
(solid curve referring to the left axis), the number of atoms 
No{t) /Nq{0) (dashed curve referring to the right axis), and the 
absolute squared overlap of the wave function with the initial 
one I J drip*{r,0)il;(r,t)\'^ (dotted curve referring to the right 
axis) with two-body dipolar loss and three-body recombina- 
tion loss. We first prepared BEG in a metastable state slightly 
below the critical point, and at f = 0, wo increased \a\ so that 
No\a\/do exceeded its critical value. The inset enlarges the 
view of the intermittent implosion. 
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FIG. 2. (a) Time evolution of the wave function \tp{r,t)\. 
A condensate of one million ^^Na atoms was prepared, and 
at t = 0, the s-wave scattering length was changed from 2.75 
nm to —1 nm. The trap frequency is a»o = 100 x 27r s~^, 
and the loss-rate coefficients are described in the text. The 
wave function at f = 0.79 is multiplied by 0.1. The inset 
shows the gray-scale image of the column density a.t t = 0.79. 
(b) Time evolution of the central height of the wave function 
\ip{r = 0,t)\ (solid curve referring to the left axis) and the 
number of atoms (dashed curve referring to the right axis). 
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